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The recent measurements of Cosmic Microwave Background temperature and polarization 
anisotropies made by the Planck satellite have provided impressive confirmation of the ACDM 
cosmological model. However interesting hints of slight deviations from ACDM have been found, 
including a 95% c.l. preference for a ’’modified gravity” structure formation scenario. In this paper 
we confirm the preference for a modified gravity scenario from Planck 2015 data, find that modified 
gravity solves the so-called Auns anomaly in the CMB angular spectrum, and constrains the am¬ 
plitude of matter density fluctuations to erg = O.SlS/Q p^g, in better agreement with weak lensing 
constraints. Moreover, we find a lower value for the reionization optical depth of r = 0.059 ± 0.020 
(to be compared with the value of r = 0.079 ± 0.017 obtained in the standard scenario), more 
consistent with recent optical and UV data. We check the stability of this result by considering 
possible degeneracies with other parameters, including the neutrino effective number, the running 
of the spectral index and the amount of primordial helium. The indication for modified gravity is 
still present at about 95% c.l., and could become more significant if lower values of r were to be 
further confirmed by future cosmological and astrophysical data. When the CMB lensing likelihood 
is included in the analysis the statistical significance for MG simply vanishes, indicating also the 
possibility of a systematic effect for this MG signal. 

PACS numbers: 98.80.-k 95.85.Sz, 98.70.Vc, 98.80.Cq 


I. INTRODUCTION 

The recent measurements of Cosmic Microwave Back¬ 
ground (CMB) anisotropies by the Planck satellite ex¬ 
periment mm have fully confirmed, once again, the ex¬ 
pectations of the standard cosmological model based on 
cold dark matter, inflation and a cosmological constant. 

While the agreement is certainly impressive, some 
hints for deviations from the standard scenario have 
emerged that certainly deserve further investigation. In 
particular, an interesting hint for ’’modified gravity” 
(MG hereafter), i.e. a deviation of the growth of density 
perturbations from that expected under General Rela¬ 
tivity (GR hereafter), has been reported in [3] using a 
phenomenological parametrization to characterize non¬ 
standard metric perturbations. 

In past years, several authors (see e.g. |3HIH]) have 
constrained possible deviations of the evolution of pertur¬ 
bations with respect to the ACDM model, by parametriz¬ 
ing the gravitational potentials $ and and their linear 
combinations. Considering the parameter E, that mod- 
ihes the lensing/Weyl potential given by the sum of the 
Newtonian and curvature potentials -|- $, the analysis 
of [5] reported the current value of Eg — 1 = 0.28 ± 0.15 
at 68% from Planck CMB temperature data, i.e. a de¬ 
viation from the expected GR null value at about two 
standard deviations. The discrepancy with GR increases 
when weak lensing data is included, bringing the con¬ 
strained value to Eg — 1 = (again, see [3]). 

This result is clearly interesting and should be further 
investigated. Small systematics may still certainly be 


present in the data and a further analysis, expected by 
2016, from the Planck collaboration could solve the issue. 
In the meantime, it is certainly timely to independently 
reproduce the result presented in [3] and to investigate 
its robustness, especially in view of other anomalies and 
tensions currently present in cosmological data. 

Indeed, another anomaly seems to be suggested by the 
Planck data, i.e. the amplitude of gravitational lensing 
in the angular spectra. This quantity, parametrized by 
the lensing amplitude Aiens as firstly introduced in m, 
is also larger than expected at the level of two standard 
deviations. The Planck-l-LowP analysis of [ 2 ] reports the 
value of Aiens = 1-22 ± O.IO at 68% c.l.. This anomaly 
persists even when considering a significantly extended 
parameter space as shown in |20j . It is therefore manda¬ 
tory to check if this deviation is in some way connected 
with the ”Eg” anomaly performing an analysis by vary¬ 
ing both parameters at the same time. This has been 
suggested but not actually done in [3]. 

Moreover, some mild tension seems also to be present 
between the large angular scale Planck LFI polarization 
data (that, alone, provides a constraint on the optical 
depth T — 0.067 ± 0.023 [ 2 ]) and the Planck HFI small- 
scale temperature and polarization data that, when com¬ 
bined with large-scale LFI polarization, shifts the con¬ 
straint to T = 0.079 ± 0.017 [ 2 ]. Since the Planck con¬ 
straints on T are mo del-dependent, is meaningful to check 
if the assumption of MG could, at least partially, resolve 
the ”t” tension. 

Another tension concerns the amplitude of the r.m.s. 
density fluctuations on scales of 8 Mpe h~^, the so-called 
erg parameter. The constraints on erg derived by the 
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Planck data under the assumption of GR and A-CDM 
are in tension with the same quantity observed by low 
redshift surveys based on clusters counts, lensing and 
redshift-space distortions (see e.g. [H] and [5]). This ten¬ 
sion appears most dramatic when considering the weak 
lensing measurements provided by the CFHTLenS sur¬ 
vey (see discussion in 0 ), which prefer lower values of 
tTg with respect to those obtained by Planck. Several so¬ 
lutions to this mild tension have been proposed, includ¬ 
ing dynamical dark energy |22j . decaying dark matter 
P5l , ultralight axions [IS] , and voids [IS] . It is there¬ 
fore timely to further check if the ’’erg tension” could be 
reconciled by assuming MG. This approach has already 
been suggested, for example, by nil. 

Finally, there are also extra parameters such as the 
running of the spectral index dns/dink, the neutrino ef¬ 
fective number N^ff (see e.g. [57]), and the helium abun¬ 
dance Yp (see e.g. [28]) that could be varied and that 
could in principle be correlated with MG. Since the val¬ 
ues of these parameters derived under A-GDM (see 0) 
are consistent with standard expectations, it is crucial 
to investigate whether the inclusion of MG could change 
these conclusions. 

This paper is organized as follows: in the next section 
we describe the MG parametrization that we consider, 
while in Section III we describe the data analysis method 
adopted. In Section IV, we present our results and in 
Section V we derive our conclusions. 


II. PERTURBATION EQUATIONS 

Let us briefly explain here how MG is implemented in 
our analysis, discussing the relevant equations. Assum¬ 
ing a flat universe, we can write the line element of the 
Friedmann-Lemaitre-Robertson-Walker (FLRW) metric 
in the conformal Newtonian gauge as: 

ds^ = a(T)^[—(1 -f 2'I')dT^ + (1 ~ 2^)dx^dxi ], (I) 

where a is the scale factor, r is the conformal time, 'k 
is the Newton’s gravitational potential, and <I> the space 
curvature 

Given the line element of the Eq. [^ we can use a phe¬ 
nomenological parametrization of the gravitational po¬ 
tentials 'k and <k and their combinations. We consider 
the parametrization used in the publicly available code 
MGCAMB [ST] [35], introducing the scale-dependent function 
fi{k,a), that modifies the Poisson equation for 'k: 

fc^Vk = —47rGa^^(fc, a)pA , (3) 


where p is the dark matter energy density, A is the co¬ 
moving density perturbation. Furthermore one can con¬ 
sider the function r]{k,a), that takes into account the 
presence of a non-zero anisotropic stress: 

’7(fc,a) = §- (4) 

We can then easily introduce the function S(fc, a), which 
modifies the lensing/Weyl potential $-|-'k in the following 
way: 

— A:^($-I-4>) = 87rGa^E(fc, a)pA , (5) 


and that can be obtained directly from p,{k, a) and r]{k, a) 
as 

Y = ^(l + V). ( 6 ) 

Of course, if we have GR then p = r/ = S = 1. 

It is now useful to give an expression for p and rj. Fol¬ 
lowing Ref. [3], we parametrize p and rj as: 


p{k,a) 

p{k,a) 


1 + /i(a) 
1 + / 2 (a) 


l + Ci{\H/kf _ 
l+\\H/kY ’ 

l + C2{\H/kf 

l + {XH/kY ’ 


(7) 

( 8 ) 


where H = d/a is the Hubble parameter, ci and C 2 are 
constants and the /i(a) are functions of time that char¬ 
acterize the amplitude of the deviation from GR. 

Again, following [3] we choose a time dependence for 
these functions related to the dark energy density: 


/i(a) = £'iir2DE(a), (9) 

where Ea are, again, constants and Ode (a) is the dark 
energy density parameter. As discussed in Ref. [5], the 
inclusion of scale dependence does not change signifi¬ 
cantly the results, we can therefore consider the scale 
independent parametrization, in which ci = C 2 = 1. 
In other words, we modify the publicly available code 
MGCAMB [3TJI35], by substituting to the original p and ry, 
the following parametrizations: 


p{k,a) = I-I-EiiQDE(a); (10) 

r]{k,a) = I-I-E 22 f^DE(a). (II) 

A detection of Ea ^ 0 could therefore indicate a depar¬ 
ture of the evolution of density perturbations from GR. 
In order to further simplify the problem, we assume a 
cosmological constant for the background evolution. 


^ In the synchronous gauge, that is the one adopted in boltzmann 
codes as CAMB m, we have: 

ds^ = a(r)^[—+ {Sij + hij)dx'^dx^] , (2) 

where hij are defined as in m- 


III. METHOD 

We consider flat priors listed in Table [I] on all the pa¬ 
rameters that we are constraining. They are: the six 
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Parameter 

Prior 

Qbh^ 

[0.005,0.1] 

Q,ch^ 

[0.001,0.99] 

0s 

[0.5,10] 

T 

[0.01,0.8] 

Us 

[0.8,1.2] 

log[10i°A,] 

[2,4] 

Ell 

[-1,3] 

E 22 

[-1.4,5] 

dris 

dink 

[-1,1] 

Neff 

[0.05,10] 

■^lens 

[0,10] 

Yp 

[0.1,0.5] 


TABLE I: External flat priors on the cosmological parameters 
assumed in this paper. 


parameters of the ACDM model, i.e. the Hubble con¬ 
stant iLg, the baryon and cold dark matter Vlch? 

energy densities, the primordial amplitude and spectral 
index of scalar perturbations, and Ug respectively, (at 
pivot scale fcg = 0.05hMpc~^), and the reionization op¬ 
tical depth r; the constant parameters of MG, En and 
E 22 ', the several extensions to ACDM model. In partic¬ 
ular we vary the neutrino effective number A^eff (see e.g. 
[27]), the running of the scalar spectral index dns/dlnk, 
the primordial Helium abundance Yp and the lensing am¬ 
plitude in the angular power spectra Hjens- We also vary 
foreground parameters following the same method of [33] 
and [ 2 ]- 

We constrain these cosmological parameters by using 
recent cosmological datasets. First of all, we consider the 
full Planck 2015 release on temperature and polarization 
CMB angular power spectra, including the large angular 
scale temperature and polarization measurement by the 
Planck LFI experiment and the small-scale temperature 
and polarization spectra by Planck HFI. We refer to the 
Planck HFI small angular scale temperature data plus 
large angular scale Planck LFI temperature and polar¬ 
ization data as Planck TT, while when we include small 
angular scale polarization from Planck HFI as Planck pol 
(see [33]). We also use information on CMB lensing from 
Planck trispectrum data (see m) and we refer to this 
dataset as lensing. Finally, we consider the weak lensing 
galaxy data from the CFHTlenS [33] survey with the pri¬ 
ors and conservative cuts to the data as described in [3] 
and we refer to this dataset as WL. 

To perform the analysis, we use our modified version, 
according to the Eqs. IT of the publicly available code 
MGCAMB [ST] [32] that modihes the original publicly code 
GAME [33] implementing the pair of functions /i(a, k) and 
77 ( 0 , fc), as dehned in [^. This code has been developed 
and tested in a completely independent way to the one 
used in [5]. 

We integrate MGCAMB in the latest July 2015 version of 
the publicly available Monte Carlo Markov Chain pack¬ 
age cosmomc |36j with a convergence diagnostic based on 
the Gelman and Rubin statistic. This version includes 
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FIG. 1: Constraints at 68% and 95% confidence levels on the 
Eg — 1 vs r plane (top panel) and on the Eg — 1 vs Hq plane 
(bottom panel) from the Planck TT and Planck pol datasets. 
The 6 parameters of the ACDM model are varied. Notice 
that Eg is different from one (dashed vertical line) at about 
95 % confidence level. A small degeneracy is present between 
Eg and r: smaller optical depths are more compatible with 
the data if Eg is larger than one (see top panel). Another de¬ 
generacy is present with the Hubble constant: larger values of 
the Hubble constant are more compatible with the considered 
data in case of Eg different from one (bottom panel). 
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FIG. 2: Constraints at 68% and 95% confidence levels on 
the Eg — 1 vs Aiens plane from the Planck TT and Planck 
pol datasets. A strong degeneracy is present between Eg and 
Aiens: larger values of Aiens are more compatible with the 
data if Eg is smaller than one. 
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Planck TT 

Planck TT ± WL Planck TT ± lensing 

Planck pol 

Planck pol ± WL Planck pol ± lensing 

Bii 

0.08l°;?l 

n 1 q+019 

0.08l°:3^ 

0.06l°;33 

—0 

O.OSlg;^® 

E 22 


1 Q+1-4 
-*-•^- 1.0 

0.4t'(® 

o.9ti;^5 

1 ^+1.3 
' - 1.0 

0At°:l 

Mo — 1 

O.Ohlg'sg 

—0 

0.05«jt 

0.04l°;23 

—0 1 5+013 

0.05l°j4 

Mo - 1 


bslJ:? 

0.31«:« 

0.6t?;g 

l.ZU_o 68 

0.261°;®® 

Eo — 1 

0.28 ±0.15 

0.34l°;l« 


0.23 ±0.13 

0.27 ±0.13 

O-lOl^i 


0.02251 ± 0.00027 

0.02263 ± 0.00026 

0.02238 ± 0.00024 

0.02237 ±0.00017 0.02243 ± 0.00017 

0.02233 ±0.00016 

Q.ch? 

0.1175 ±0.0024 

0.1159 ±0.0022 

0.1171 ±0.0021 

0.1188 ±0.0016 

0.1180 ±0.0015 

0.1185 ±0.0014 

Ho 

68.5 ± 1.1 

69.2 ± 1.1 

68.47 ±0.99 

67.78 ± 0.71 

68.15 ±0.69 

67.83 ±0.66 

T 

0.065 ±0.021 

n ofii 

U.UOi_o.o23 

0.050 ±0.019 

0.059 ±0.020 

0.054 ±0.019 

0.045 ±0.017 

ris 

0.9712 ±0.0071 

0.9754 ±0.0067 

0.9706 ± 0.0062 

0.9668 ±0.0051 

0.9689 ± 0.0050 

0.9668 ± 0.0047 

eg 

U.01U_Q Q52 

0.787lHf9 

0.802t°°fg 

n 81 5+0 032 

U.OlO_Q Q43 

0-788in35 

0.803 ±0.031 


TABLE II: Constraints at 68% c.l. on the cosmological parameters assuming modified gravity (parametrized by -Bn and B 22 ) 
and varying the 6 parameters of the standard ACDM model. 



FIG. 3: Constraints at 68% and 95% conhdence levels on 
the So — 1 vs ALff plane from the Planck TT and Planck 
polarization datasets. Notice that So is different from unity 
(dashed vertical line) at about the 95 % confidence level. A 
small direction of degeneracy is present between So and N^s'- 
larger Nsa are more compatible with the data if So is larger 
than one in case of the Planck TT dataset. 


the support for the Planck data release 2015 Likelihood 
Code |33j (see http://cosmologist.info/cosmomc/) 
and implements an efficient sampling using the fast/slow 
parameter decorrelations ISZj. 


IV. RESULTS 

We first report the results assuming a modified grav¬ 
ity scenario parametrized by rj and ^ and varying only 
the 6 parameters of the standard ACDM model. The 
constraints on the several parameters are reported in Ta¬ 
ble 1^ When comparing the first and second column of 
our table, we see a complete agreement with the results 
presented in the first and third column of Table 6 of [3] . 
Namely we find evidence at ~ 95% c.l. for different 
from zero for the Planck TT dataset, and this indication 
is further confirmed when the WL dataset is included. 

As fully discussed in [33], the Planck polarization HFI 
data at small angular scales fails to satisfy some of the 
internal checks in the data analysis pipeline. The re¬ 
sults obtained by the inclusion of this dataset should 
therefore be considered as preliminary. We report the 
constraints from the Planck pol dataset in columns 4-6 
in Table |TT| As we can see, the small angular scale HFI 
polarization data improves the constraints on Eg, also 
slightly shifting its value towards a better compatibility 
with standard ACDM. We can see however that the in¬ 
clusion of small angular scale polarization does not alter 
substantially the conclusions obtained when using just 
the Planck TT dataset. 

Considering just the Planck TT dataset, it is interest¬ 
ing to note that in this modified gravity scenario, the 
Hubble constant is constrained to be Hq = 68.5 ± 1.1 
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Planck TT 

Planck TT ± WL Planck TT ± lensing 

Planck pol 

Planck pol ± WL Planck pol ± lensing 

Ell 


—0 1 

o.o8l°:f3 

0.07l°;33 


0.06l°;33 

E 22 

0.81- 

-*-•^-1.3 

o.8ll:° 

0.7l);° 

1.4± 1.2 

o.8l);J 

Mo — 1 

0.04l°;23 


0-06lgji 

0.05l°i3 

—0 

^•-L^_0.33 

0-04l°j2 

Mo - 1 

0 . 61 ?;^ 

1 n+1.0 

-‘-•'^-0.9 

o.5l?:I 


0.95 ±0.81 


Eg — 1 

U.Z1_Q 21 

U.^^_0,22 

0 . 21 «; 1 ® 

0.19lo.;8 

^•^^-0.18 

U.zz,_o 16 


0.02259 ± 0.00029 

0.02273 ± 0.00028 

0.02231 ±0.00026 

0.02239 ± 0.00017 0.02246 ± 0.00017 

0.02229 ±0.00016 

Q.ch? 

0.1169 ±0.0025 

0.1152 ±0.0023 

0.1180 ±0.0025 

0.1187 ±0.0016 

0.1177 ±0.0015 

0.1191 ±0.0015 

Ho 

68.8 ± 1.2 

69.6 ± 1.1 

68.1 ± 1.2 

67.82 ± 0.73 

68.26 ±0.69 

67.59 ±0.70 

T 

0.059l°:°23 

0.054 ±0.021 

0.059 ±0.021 

0.056 ±0.020 

0.0491^22 

0.057 ±0.021 

ris 

0.9730 ± 0.0073 

0.9772 ± 0.0068 

0.9687 ±0.0070 

0.9671 ±0.0050 

0.9694 ± 0.0049 

0.9656 ± 0.0049 

eg 

0 «n7+0 033 
U.OU/ _0.049 

n 7 Q 9 + 0.025 

u. 1 oz,_Q Q 3 g 

0.813l°:°f6 

0 81 

U.O-LO_q Q 44 

n 7Qfi;+0-023 

U.78d_q Q 35 

0 81 A+O '^^i 
U-Ol^_0 046 

-^lens 

i.uy_Q 13 

1 iq + 0.10 
-L--LO_o.i4 

0.924l°°fg 

1 Q^+O.OS 

1 Q7+0.09 
- 0.11 

U.yi^_0_078 


TABLE III: Constraints at 68 % c.l. on the cosmological parameters assuming modified gravity (parametrized by En and E 22 ) 
and varying the 6 parameters of the standard ACDM model plus Aiens- 


at 68 % C.L, i.e. a value significantly larger than the 
Hq = 67.3 ± 0.96 at 68 % c.l. reported by the Planck 
collaboration assuming ACDM. Combining the Planck 
TT dataset with the HST prior of Hq = 73.0 ± 2.4 from 
the revised analysis of [32] as in [ 33 ] we found indeed an 
increased evidence for MG, with Eq — 1 = 0.33^° at 
68 % C.L. 

Moreover, the amplitude of the r.m.s. mass density 
fluctuations CTg in our modified gravity scenario is con¬ 
strained to be cTg = 0 . 816 lg;Q 52 at 68 % c.l., i.e. a value 
significantly weaker (and shifted towards smaller values) 
than the value of erg = 0.829 ± 0.014 at 68 % c.l. re¬ 
ported by the Planck collaboration again under ACDM 
assumption. 

Considering the Planck pol dataset, the value of the op¬ 
tical depth is also significantly smaller in the MG scenario 
(r = 0.059 ± 0.020 at 68 % c.l.) respect to the value ob¬ 
tained under standard ACDM model of t = 0.078±0.019 
at 68 % C.L, i.e. reducing the tension with the Planck LFI 
large angular scale polarization constraint. Interestingly, 
a smaller value for the optical depth of r ~ 0.05 is in 
better agreement with recent optical and UV astrophys- 
ical data (see e.g. |44H46] 1 and the reionization scenarios 
presented in |38|. A value of r > 0.07 could imply un¬ 
expected properties for high-redshift galaxies. Assuming 
an external gaussian prior of r = 0.05 ± 0.01 (at 68 % 


C.L.) as in |48j that would consider in a conservative way 
reionization scenarios where the star formation rate den¬ 
sity rapidly declines after redshift 2 ; ~ 8 as suggested by 
m, we find that the Planck TT dataset provides the 
constraint Eg — 1 = 0.30 ± 0.14 at 68 % c.l., i.e. further 
improving current hints of MG. In this respect, future, 
improved, constraints on the value of r from large-scale 
polarization measurements as expected from the Planck 
HFI experiment will obviously provide valuable informa¬ 
tion. 

The degeneracies between Eg, Hq and r can be clearly 
seen in Figure [T] where we show the constraints at 68% 
and 95% confidence levels on the Eg — 1 vs r plane (top 
panel) and on the Eg — 1 vs Hq plane (bottom panel) 
from the Planck TT and Planck pol datasets. As we 
can see, a degeneracy is present between Eg — 1 and r: 
smaller optical depths are more compatible with the data 
if Eg is larger than one (see top panel). As discussed, a 
second degeneracy is present with the Hubble constant: 
larger values of the Hubble constant are more compatible 
with the considered data in case of Eg different from one 
(Bottom Panel). 

As already noticed in |3] and as we will discuss in the 
next paragraph, the indication for MG from the Planck 
data is strictly connected with the Aie„s anomaly, i.e. 
with the fact that Planck angular spectra show ’’more 
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Planck TT 

Planck TT ± WL Planck TT ± lensing 

Planck pol 

Planck pol ± WL Planck pol ± lensing 
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n ny+O.Sl 

U.U ( _o,73 

—0 1 ‘^+'^■ 2 '^ 
U.iO_o. 5 g 
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0-0711®^ 

-0-2lll?^ 

0-08111? 

E 22 
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2 -ill;g 

0.511° 

0-9lli 

1-7511:^ 

0.4l?;® 

Mo — 1 

0 - 05 l°j 2 

-0.09l°i| 

O.OOll.fs 

0.05li;® 

—0 1 

U.iO_Q 33 

0.06111? 

770 - 1 

0.96 ± 1.1 

i-5lS;| 

0-311? 

0.59l?;^ 

-1 99+0.96 
-L-^^_0.69 

0-24ll®« 

Eo — 1 

0.36 ±0.18 

0.45l°l) 

U.lZ_o_i 4 

U.ZO_Q 25 

0.2811,;? 

o-ioll?? 


0.022941^0^063 

n nooQQ+0.00048 

U.UZOZO_Q QQQg3 

n 

U.UZZOZ_Q QQQ43 

0.02234 ± 0.00025 

0.02244 ± 0.00026 

0.02224 ±0.00024 


0.1202 ±0.0041 

U.1Z1U_Q QQ4g 

0.1185 ±0.0039 

0.1184 ±0.0030 

0.1181 ±0.0030 

0.1173 ±0.0030 

Ho 

72.011® 

74.711® 

69.711® 

67.6 ± 1.6 

68.2i;° 

67.1 ±1.5 

T 

U.UI ^_o ,026 

0.073 ±0.024 

0.052ll°"2® 

0.05911°;? 

n Qt::o+0.019 
U.UJO_Q Q21 

0-04411?;? 

Us 

U.yyu_Q Q 25 

l.UU^_Q Q 25 

0.97711.°;? 

0.9655 ± 0.0097 

0.969 ± 0.010 

0.9625 ± 0.0092 

0-8 

n qo' 7 + 0-033 
u.oz i _ 0.062 

n SI n+O-028 
U.»i^_ 0.054 

0.808iroi| 

n SI 4 +0 032 
(J.8i4_Q Q4g 

n ' 7 QO+ 0.022 

U. I oo_Q 033 

0.7991?;??; 

Aeff 

0 ^1 +0.36 

Q ^0+0.35 

0.00_o 43 

3.19ll.“ 

3.02 ±0.20 

3.06 ± 0.21 

2.95 ±0.19 


TABLE IV: Constraints at 68 % c.l. on the cosmological parameters assuming modified gravity (parametrized by En and E 22 ) 
and varying the 6 parameters of the standard ACDM model plus ALn. 


lensing” than expected in the standard scenario. It is 
therefore not a surprise that when the Planck lensing 
data (obtained from a trispectrum analysis) that is on the 
contrary fully compatible with the standard expectations 
is included in the analysis the indication for modified 
gravity is significantly reduced to less than one standard 
deviation, as we can see from the third column of Table 
[nj On the other hand, when weak lensing data from the 
WL dataset is included, the indication for MG increases, 
with So ~ 1 larger than zero at more than 95% c.l.. 


In Tables pHl |IVI [Vl [VI| we report constraints assuming 
one single parameter extension to ACDM. In particular, 
we report constraints when adding as an extra parameter 
the lensing amplitude Ajens (Table [TTl| , the neutrino ef¬ 
fective number Ngg (Table [Iv| ), the running of the scalar 
spectral index dns/dlnk (Table [v| and, finally, the He¬ 
lium abundance Yp (Table [Vl|. 


As expected, there is a main degeneracy between the 
Aiens parameter and Sq, as we can clearly see in Figure 
where we report the 2D posteriors at 68% and 95% 
c.l. in the Sq — 1 vs Aiens plane from the Planck TT 
and Planck pol datasets. In practice, the main effect of 
a modified gravity model is to enhance the lensing signal 
in the angular power spectrum. The same effect can be 
obtained by increasing Aiens and some form of degeneracy 
is clearly expected between the two parameters. As we 


see from the results in Table |III[ the value of the Aiens 
parameter, when MG is considered, is Ajens = 1-091°'ijg, 
fully consistent with I, while for the standard ACDM 
the constraint is Aiens = When 

also varying Aiens we found that the Planck pol datasets 
constraint the optical depth to t = 0.056 ± 0.020 at 68% 
c.l. 


On the other hand, by looking at the results in Ta¬ 
bles W V VI we do not see a significant degeneracy 
between the MG parameters and the new extra param¬ 
eters. A small degeneracy is however present between 
So and the effective neutrino number Nes- We see from 
Table IV that Planck TT data provides the constraint 
iVefi = 3.4llQ;4g at 68% c.l. that should be compared 
with A"eff = 3.131 q'34 at 68% c.l. from the same dataset 
but assuming the standard ACDM model. While the pos¬ 
sibility of an unknown ’’dark radiation” component (i.e. 
iVeff > 3.046, see e.g. pHIIiT)] ! is therefore more viable in 
a MG scenario, it is however important to note that when 
adding polarization data the constraint on the neutrino 
number is perfectly compatible with the expectations of 
the standard three neutrino framework. The constraints 
at 68% and 95% c.l. in the Sq — 1 vs Ngg planes are 
reported in Figure 

We also consider the possibility of a running of the 
scalar spectral index dns/dlnk. Results are reported in 
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TABLE V: Constraints at 68% c.l. on the cosmological parameters assuming modified gravity (parametrized by -Eii and E 22 ) 
and varying the 6 parameters of the standard ACDM model plus dns/dlnk. 


Table [V| and we find no degeneracy with MG parameters. 
The Planck TT constraint of dns/dlnk = —O.OOTSto^oose 
at 68% c.l. is almost identical to the value dns/dlnk = 
—0.0084 ± 0.0082 at 68% c.l. obtained using the same 
dataset but assuming standard ACDM. 

We also considered variations in the primordial he¬ 
lium abundance Yp since it affects small angular scale 
anisotropies. Our results are in Table |VI[ The Planck 
TT constraint is found to be Yp = 0.258 ± 0.023 at 68% 
C.L, slightly larger than the standard ACDM value of 
Yp = 0.252 ± 0.021 at 68% c.l. obtained using the same 
dataset. While a larger helium abundance is in better 
agreement with recent primordial helium measurements 
of [H], it is important to stress that the inclusion of po¬ 
larization yields a constraint that is almost identical to 
the one obtained under ACDM. The constraints at 68% 
and 95% c.l. in the Eg ~ 1 vs dns/dlnk and Eq — 1 vs 
Yp planes are reported in Figure]^ 


V. CONCLUSIONS 

In this paper, we have further investigated the current 
hints for a ’’modified gravity” scenario from the recent 
Planck 2015 data release. We have confirmed that the 
statistical evidence for these hints, assuming the conser¬ 


vative dataset of Planck TT, is, at most, at ^ 95% c.l., 
i.e. not extremely significant. The statistical significance 
increases when combining the Planck datasets with the 
WL cosmic shear dataset. Indeed, the Planck dataset 
seems to provide lower values for the erg parameter with 
respect to those derived under the assumption of GR and 
A-CDM. 

If future astrophysical or cosmological measurements 
will point towards a lower value of the optical depth of 
T ~ 0.05 or of the r.m.s. amplitude of mass fluctuations 
of (Tg ^ 0.78 then the current hints for modified gravity 
could be further strenghtened. 

However it also important to stress that when the CMB 
lensing likelihood is included in the analysis the statisti¬ 
cal significance for MG simply vanishes. 

We also investigated possible degeneracies with extra, 
non-standard parameters as the neutrino effective num¬ 
ber, the running of the scalar spectral index and the pri¬ 
mordial helium abundance showing that the results on 
these parameters assuming ACDM are slightly changed 
when considering the Planck TT dataset. Namely, under 
modified gravity we have larger values for the neutrino 
effective number, N^ff = at 68% c.l., and for 

the helium abundance, Yp = 0.258 ± 0.023. at 68% c.l.. 
However, the constraints on these parameters are practi¬ 
cally identical those obtained under GR when including 
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TABLE VI: Constraints at 68 % c.l. on the cosmological parameters assuming modified gravity (parametrized by En and E 22 ) 
and varying the 6 parameters of the standard ACDM model plus Yp. 


the Planck HFI polarization data. 

We have clearly shown that the slight Planck hints 
of MG are strongly degenerate with the anomalous 
lensing amplitude in the Planck CMB angular spectra 
parametrized by the Ajens parameter. Indeed, the Aiens 
anomaly disappears when MG is considered. Glearly, un¬ 
detected small experimental systematics could be the ori¬ 
gin of this anomaly. However our conclusions are that 
modified gravity could provide a physical explanation, 
albeit exotic, for this anomaly that has been pointed out 
already in pre-Planck CMB datasets [35], was present in 
the Planck 2013 data release m and seems still to be 
alive in the recent Planck 2015 release [2| 

An extra parameter we have not investigated here is 
the neutrino absolute mass scale Sm^. Since MG is de¬ 
generate with the Aiens we expect that in a MG scenario 
current constraints on the neutrino mass from CMB an¬ 
gular power spectra should be weaker. However a more 
detailed computation is needed and we plan to investi¬ 
gate it in a future paper ([52]). 

During the submission process of our paper, another 
paper appeared [53] , claiming an indication for MG from 


cosmological data. The dataset used in that paper is 
completely independent from the one used here and the 
MG parametrization is also different. Clearly a possi¬ 
ble connection between the two results deserves future 
investigation. 
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